Stability of nonstationary states of spin-2 Bose-Einstein condensates 
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The dynamical stability of nonstationary states of homogeneous spin-2 rubidium Bose-Einstein 
condensates is studied. The states considered are such that the spin vector remains parallel to the 
magnetic field throughout the time evolution, making it possible to study the stability analytically. 
These states are shown to be stable in the absence of an external magnetic field, but they become 
unstable when a finite magnetic field is introduced. It is found that the growth rate and wavelength 
of the instabilities can be controlled by tuning the strength of the magnetic field and the size of the 
condensate. 
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I. INTRODUCTION 



The physics of F = 2 spinor Bose-Einstein condensates 
(BECs) started to gain the attention of both theorists 
and experimentalists during the last decade. The inter- 
est is motivated by the structure of F = 2 condensates: 
being more complex than that of _F = 1 condensates, 
it makes possible properties and phenomena which are 
not present in an F = 1 system. One example of this 
can been seen in the structure of the ground states. The 
energy functional of an F = 2 condensate is character- 
ized by one additional degree of freedom compared to the 
F = 1 case. This leads to a rich ground state manifold 
as now there are two free parameters parametrizing the 
ground states P, Q • This should be contrasted with an 
F = 1 condensate, where the ground state is determined 
by the sign of the spin-dependent interaction term 0] . 
Another difference can be seen in the structure of topo- 
logical defects. It has been shown that non-commuting 
vortices can exist in an F = 2 condensate ^ , while these 
are not possible in an F = 1 BEG 0,13. The topological 
defects of F = 2 condensates have been studied further 
in Refs. Experimental studies of F = 2 BECs 

have been advancing in the past ten years. Experiments 
on F = 2 ^''Rb atoms cover topics such as spin dynamics 
p^l - [l3 |. creation of skyrmions [l^, and spin-dependent 
inelastic collisions [l^. An F = 2 spinor condensate of 
^^Na atoms has been obtained experimentally Q, but it 
has a much shorter lifetime than F = 2 rubidium con- 
densates. 

In this work, we study the dynamical stability of non- 
stationary states of homogeneous F = 2 spinor conden- 
sates. The stability of stationary states has been ex- 
amined in Refs. whereas that of nonstationary 
states has received only little attention. Previous studies 
on the topic concentrate on F = 1 condensates p^| - |23j . 
Here we extend the analysis of Ref. [2^ to an F = 2 
rubidium condensate and present results concerning the 
magnetic field dependence of the excitation spectrum and 
stability. Although we concentrate on the stability of 
®^Rb condensates, many of the excitation spectra and 
stability conditions given in this article are not specific 
to rubidium condensates but have a wider applicability. 



We show that, in comparison with an F = 1 system, the 
stability analysis of an F = 2 condensate is consider- 
ably more complex. This is partly due to the presence of 
a spin-singlet term in the energy functional of the latter 
system, but the main reason for the increased complexity 
is seen to be the much larger number of states available 
in an F = 2 condensate. 

This article is organized as follows. Section |IT] intro- 
duces the system and presents the Hamiltonian and equa- 
tions of motion. In Section IIIII the Bogoliubov analysis 
of nonstationary states is introduced. This method is 
applied to study the stability both in the presence and 
absence of a magnetic field. In this section it is also de- 
scribed how Floquet theory can be used in the stability 
analysis. In Section ITVl the stability is studied under the 
(physically motivated) assumption that one of the inter- 
action coefficients vanishes. Finally, section |V] contains 
the concluding remarks. 



II. THEORY OF A SPIN-2 CONDENSATE 

The order parameter of a spin-2 Bose-Einstein con- 
densate can be written as i/j — (7/;2, V'lj V'Oj V'-ij ^-2)"^, 
where T denotes transpose. The normalization is 
Sm=-2 iV'mP = where n is the total particle den- 
sity. We assume that the trap confining the condensate 
is such that all the components of the hyperfine spin can 
be trapped simultaneously and are degenerate in the ab- 
sence of magnetic field. This can be readily achieved in 
experiments [2^ . If the system is exposed to an external 
magnetic field which is parallel to the z axis, the energy 
functional reads 



EH'] = / dr 



(1) 



where F = {Fx,Fy,Fz) is the (dimensionless) spin oper- 
ator of a spin-2 particle. & describes singlet pairs and is 
given by = 2V'2V'-2 — 2V'iV'-i +V'o- 1^ can also be writ- 
ten as O = ^'^e~'^'^^yip- The single-particle Hamiltonian 
h reads 



2m 



C/(r)-Ai-pF, +gF2 



(2) 
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Here U is the external trapping potential, ^ is the chemi- 
cal potential, and p = —g^b^B is the linear Zeeman term. 
In the latter g is the Lande hyperfine g-factor, /ib is the 
Bohr magneton, and B is the external magnetic field. 
The last term in Eq. ^ is the quadratic Zeeman term, 
q = —{g^BB)^/Ehf, where Em is the hyperfine splitting. 
The sign of q can be controlled experimentally by using 
a linearly polarized microwave field [25i] . In this article 
we consider both positive and negative values of q. 

The strength of the spin-independent interaction is 
characterized by go = 47rft^(4a2 -I- 3a4)/7m, whereas 
gi = 47rfi,^(a4 — a2)/7m and 52 = 47r?i^[(ao — 04)75 — 
2(a2 — 04)/?] describe spin-dependent scattering. Here 
Of is the s-wave scattering length for two atoms collid- 
ing with total angular momentum F. In the case of ® ''Rb, 
we calculate 170 using the scattering lengths given in Rcf. 
[l| , and 52 and 34 are calculated using the experimentally 
measured scattering length differences from Ref. [2^ . 

Two important quantities characterizing the state ip 
are the spin vector 



f(r) 



n{r) 



(3) 



and the magnetization in the direction of the magnetic 
field 



AL = 



Jdrn{r)f,{r) 
Jdrn{r) 



(4) 



The length of f is denoted by /. For rubidium the mag- 
netic dipole-dipole interaction is weak and consequently 
the magnetization is a conserved quantity. The Lagrange 
multiplier related to the conservation of magnetization 
can be included into p. The time evolution equation ob- 
tained from Eq. ([T]) is 

zh^^ij = (5) 



where 



H[i^] = h + .goV/ V + .91 (F) • F + g2Qf. (6) 



Here T 



"C* is the time-reversal operator, where 



C is the complex conjugation operator. 



III. STABILITY OF NONSTATIONARY STATES 
WHEN 32 / 



that Q vanishes at all times, and we now study this case. 
We define a state 



* 



V y 



-!</.< 2. (7) 



For this state 6 = 0, (F^) ^ (Fy) = 0, and (F.) 
fz- Furthermore, the populations of the state il'2;-i re- 
main unchanged during the time evolution determined 
by the Hamiltonian ([6|). Consequently, G = through- 
out the time evolution. The state ip2;-i with = 0, 
called the cyclic state, is a ground state at zero mag- 
netic field The creation of vortices with fractional 
winding number in states of the form 'i/'2;-i has been 
discussed in 0. The stability properties of the state 
= V^iO, V^TTz, 0, 0, VT^7r)^/\/3 are similar to 
those of V-'2:-i and will therefore not be discussed sepa- 
rately. 

The Hamiltonian giving the time evolution of 'ip2:-i is 



H[ip2--i] = gon - /-i + iginfz - p)Fz + qP. 



(8) 



where we have set f7 = as the system is assumed to 
be homogeneous. This is of the same form as the Hamil- 
tonian of an _F = 1 system discussed in [2^. The time 
evolution operator of V'2;-i is given by 



(9) 



We analyze the stability in a basis where the state 
"02;-! is time-independent. In the new basis, the energy 
of an arbitrary state (/) is given by [23l | 

= i?[?72;-10] + lh{<^\ (^?^2T-l) '^2;-10>, (10) 

and the time evolution of the components of (\) can be 
obtained from the equation 

We replace — ^ '02;-i+(5V-' in the time evolution equation 
([TT|) and expand the resulting equations to first order in 
bilj. The perturbation = {S'ipi, Sipi, Sipo, Sijj-i, 5'4'-2)'^ 
is written as 



The stability analysis is performed in a basis where the 
state in question is time-independent. This requires that 
the time evolution operator of the state is known. As 
we are interested in analytical calculations, an analytical 
expression for this operator has to be known. In order 
to calculate the time evolution operator analytically, the 
Hamiltonian has to be time-independent. In particular, 
the singlet term should not depend on time. This is 
clearly the case if the time evolution of the state is such 



where j = —2,-1,0,1,2. Straightforward calculation 
gives the differential equation for the time evolution of 
the perturbations as 



dt Ivk/ ^' I Vk 



-82.-1 = 



X -Y 
Y* -X* 



(12) 
(13) 
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where Uk = (u2;k, ui;k, wo;k, w-i;k, w-2;k)'^, Vk is defined 
similarly, and the 5x5 matrices X and Y arc 

^ = efc+50|V'2;-l>(V'2;-l|+.9lE I V'^ - 1 W) (V^- 1 W I 



-2g2|V'2:-l(i)>(^2;-lWI, 



(14) 



f=go|^2;-l)(V2Vll+5lE l^2;-lW>((^i-l)*W|- 



(15) 



The eigenvectors can be calculated straightforwardly, 
see Appendix El The eigenvectors, like the eigenvalues, 
are independent of 92 ■ The perturbations corresponding 
to the eigenvectors of can be written as 



1,2,3,4 



(r, k; t) = Ci,2,3,4(r, k; t) 'ip2-- 



(22) 



where Cj's include all position, momentum and time de- 
pendence. These change the total density of the conden- 
sate and are therefore called density modes. 



Here we have defined 



2m 

rt 



(16) 



^2;-lW = C^2;-l(*)-^jt^2;-l(i)V'2;-l, J=X,y,Z, (17) 



V'2;-lW = Ui.,it)e''^'yU2;-im2, 



(18) 



In the rest of the article we call the operator determining 
the time evolution of the perturbations the Bogoliubov 
matrix. In the present case, -62:-! is the Bogolibov ma- 
trix of V'2;-i- It is possible to write i?2:-i as a direct sum 
of three operators 

B2;-l(i) = S2V1 ® ^2Vl(^) ® ^2;'-l(0, (W) 

Bl[-i = -(Bl-ir- (20) 

Here is a time-independent 4x4 matrix and 

is a time-dependent 3x3 matrix. The bases in which 
these operators are defined are given in Appendix [A] 
The time-dependent terms of -Bf ^i are proportional to 
^±ikqt/h^ where fc = 2,4 or 6, and consequently the sys- 
tem is periodic with minimum period T = ■nh/q. Hence 
it possible to use Floquet theory to analyze the stability 
of the system [23| • In the following we first calculate the 
eigenvalues of B2._i, then those of -Sf-.i and B2 _i 
the case q = 0, and finally we discuss the general case 
q ^ using Floquet theory. 



A. Eigenvalues of B2; 



First we calculate the eigenvalues and eigenvectors of 
_B|._]^. This operator is independent of q. The eigenval- 
ues are 



huji 



2,3,4 



± 



Efc (cfe + gon + 5in(2 -I- /^) 



± n^[g„ - .91 (2 + /,)]2 + 4go5i/|)] (21) 



Here we use a labelling such that ++, — h, H — , and 

correspond to a;i,W2,W3, and W4, respectively. Now cJi.2 
have a non- vanishing imaginary part only if 170 and gi are 
both negative, while ^3^4 have an imaginary component 
if go and gi are not both positive. Consequently, these 
modes arc stable for rubidium for which go^gi > 0. 



B. Eigenvalues of B^-^-i and Bf;-! at q = 

In the absense of an external magnetic field -B|._i is 
time-independent. The eigenvalues of Sf -i can be ob- 
tained from those of -Bf .i by complex conjugating and 
changing the sign. For this reason we give only the eigen- 
values of -Bf -i- 



1 r 

2 



232", 



(23) 



ginf, ± v/(2efc - gmf.Y + l&gmek . (24) 



The eigenvalues tvjjQ and hujj have a non- vanishing com- 
plex part if 51 < 0. For rubidium all eigenvalues are 
real. There are two gapped excitations: at gfc = we 
get hu)^ = 2.92?^ and fnuQ (hcuj) = ginf^ if gifz > 
(.9i/z < 0)- The eigenvectors are given in Appendix A. 
The corresponding perturbations become 







^^^(r,k;t) = C5 







\-VTTf: 

I 



\ 



(25) 



<5V''^(r,k;t)=EC, 



6,7 



3*''-'-ginV(2 -/,)(! + /,) 

""'■""V^(efc + 25i^-^6,7) 


X e"^"-5in(2-/,), / 

(26) 

where Cs.g.r are functions of r, k and t. These modes 
change both the direction of the spin and magnetization 
and are therefore called spin-magnetization modes. 



C. Non- vanishing magnetic field 

If 9 7^ 0, the stability can be anal yze d using Floquet 
theory due to the periodicity of i?2 _i Hi- We denote the 
time evolution operator determined by B2._i by J/l-i- 
According to the Floquet theorem (see, e.g., [131), ^^^l:-! 
can be written as 



-UK 



(27) 
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where M is a periodic matrix with minimum period T 
and Af (0) = I, and K is some time-independent ma- 
trix. At times t ~ nT, where n is an integer, we get 
f/|._i(7iT) = e"™^-^. The eigenvalues of K determine 
the stabihty of the system. We say that the system is 
unstable if at least one of the eigenvalues of K has a pos- 
itive imaginary part. We calculate the eigenvalues {fko} 
of K from the equation 



In A 



(28) 



where {A} are the eigenvalues of U2.^i{T). 




q=-l 



q=-5 




FIG. 1. The positive imaginary part uf related to U2--i{T) for 
different values of the magnetic field parameter q (measured 
in units of |(;i|7i). The units of et and a;' are \gi\n and \gi\n/h, 
respectively. Note that the scales of the and cu' axes are 
not equal in the top and bottom rows. 

We plot w' for several values of the magnetic field in 
Fig. [ij By comparing to the case of a rubidium con- 
densate with 92 =0, we found that the instabilities are 
essentially determined by gi, the effect of g2 is negligible. 
The eigenvectors of U2--i{T) correspond to perturbations 
which affect both spin direction and magnetization. With 
the help of numerical results we find that a good fitting 
formula is given by 



(29) 



.lm{^(efc -2q + g^nf - -|(/. - 2)(/, + l)gin|}, 
g > 0. (30) 



We see that for q> Q the fastest-growing instability is lo- 
cated approximately at efc = max{0, 2q — gin\ regardless 
of the value of . For q <Q the location of this instability 
becomes magnetization dependent and is approximately 
given by eu = max{0, |g| — 5(2 — fz)\gi\n/&}. The values 
of Efe corresponding to unstable wavelengths are bounded 
above approximately by the inequality efc < (Sjgl -I- q)/2. 
Therefore, the state ■02: -i is stable if the condensate is 
smaller than the shortest unstable wavelength 



A2:-l 



2i:h 



vMSM + g) 



(31) 



At q = the system is stable regardless of its size. 

Fig. [1] shows that the shape of the unstable region 
depends strongly on the sign of q. This can be understood 
qualitatively with the help of the energy functional of Eq. 
([T]). We choose | — 1) to be the initial state of the system 
and assume that the condensate breaks into domains of 
length L such that each domain is either in state |0) or 
in state | — 2) and that neighboring domains are in a 
different state. Then the energy of the initial state is 
gin/2 + q (dropping constant terms), while the energy of 
the final configuration is gin + 2q. If gi,q > 0, domain 
formation is suppressed for energetic reasons. If, on the 
other hand, gi > and q < 0, this process is enhanced 
when q < —gin/2. This leads to an instability near = 
— 1 if g is negative, but makes the state stable in this 
region if q is positive. 



IV. STABILITY OF NONSTATIONARY STATES 
WHEN ff2 = 

For rubidium the value of g2 is small in comparison 
with go and gi. Consequently, it can be assumed that 
this term has only a minor effect on the stability of the 
system. This assumption is supported by the results of 
the previous section. In the following we will therefore 
study the stability in the limit 52 = 0. This makes it 
possible to obtain an analytical expression for the time 
evolution operator also for states other than V'2:-i- First 
we discuss a state that has three nonzero components, 
and then two states that have two nonzero components. 



A. Nonzero tp2, i^o, and tp-2 
We consider a state of the form 



/V2 - 2po + h\ 
' ' 



V'2;0;-2 = — 



26'"^ 



VV2 - 2po - /;/ 



IM <2-2po (32) 



For this (F^) = (Fy) = and (F^) = nf^. The Hamilto- 
nian and time evolution operator of this state are given 



5 



by Eqs. ^ and ([9]), respectively. The equations de- 
termining the time evolution of the perturbations can be 
obtained from Eqs. (fT4|) and (fT5|) by replacing V'2;-i with 
i^2;0;-2 and setting g2 = 0. In this way, we obtain a time 
dependent Bogoliubov matrix -B2;0;-2 which is a func- 
tion of the population of the zero component, po- The 
Bogoliubov matrix can now be written as 



^2;0;-2(i) 



'^2;0;-2 



■82;0;-2(0: 



(33) 



where i32 o -2 time independent and i?|.o _2 periodic 
in time with period T = Tr/q. The bases in which these 
operators are defined are given in Appendix B. The eigen- 
values of S2 -2 ^"^^ 



If g < 0, the state is stable regardless of the size of the 
condensate. In Fig. [5]we plot the positive imaginary part 
of the eigenvalues ((36|) for various values of q. 



Po=0 



Po=0.l 



Po=0.5 




2 4 6 8 10 



2 4 6 8 10 



2 4 6 8 10 



huj 



3,4,5,6 



Sfc + efe[.9o + 4.9i(l-po)] 



± <^kn^y[go - 4gi(l - po)]^ + 4go.9i/| 



(34) 



1/2 



(35) 



Here — +, H — , and correspond to wa, W4, ws, and 

We, respectively. These eigenvalues are always real if go 
and gi are positive. From the eigenvectors given in Ap- 
pendix [B] we see that ^3^4 are density modes and aji.2 
and uj^^Q are magnetization modes. All these are gapless 
excitations. Note that the eigenvalues are independent 
of 61. 

We discuss next the stability properties determined by 
We consider first the special case po = and 
proceed then to the case po > 0. 



r4 

-°2;0;-2- 



1. Stability at po ~ 

In the case po = a complete analytical solution of 
the excitation spectrum can be obtained. In Appendix 
IB] we show that by a suitable choice of basis the time 
dependence of the Bogoliubov matrix can be eliminated. 
The eigenvalues are 



1 

.10-2 



± ginfz + 6q 



± \/4(efc + gm - 3<7)2 - (4 - f^){giny 



(36) 



These are gapped excitations and correspond to spin- 
magnetization modes (see Appendix [B|) . If gi > 0, these 
eigenvalues have a non-vanishing complex part when 
3g— 2gin < efe < 3g. This is possible only if q is positive. 
The location of the fastest-growing unstable mode, de- 
termined by efc = max{0, 3q — gin}, is independent of fz- 
The maximal width of the unstable region in the efe di- 
rection, obtained at fz = 0, is 2|(7i|n. The state is stable 
if it is smaller than the size given by 



A 



2:0: 



2TTh 

y/Gmq ' 



q>0. 



(37) 



FIG. 2. The positive imaginary part tj' related to the eigenval- 
ues (|36|l and to U2;0;-2{T) for various values of the quadratic 
Zeeman term q and population po. The unit of q and is 
\gi\n, and that of is \gi\n/h. We have chosen = as this 
choice gives the fastest-growing instabilities and the smallest 
size of a stable condensate. In the top row the dashed, dotted, 
and solid lines correspond to q = 1, 2, 3, respectively, while in 
the bottom row they correspond to q = —1, —2, —3, respec- 
tively. We have set 6 — in i/)2;0;-2 as the stability was found 
to be independent of 6. 



2. Stability when po > 

In the case po > the stability can be studied using 
Floquct theory. The stability properties can be shown 
to be independent of the sign of fz- At g = the op- 
erator i3|.o _2 is time- independent. The eigenvalues can 
be obtained analytically but are not given here. The 
eigenvalues show that in the absence of magnetic field 
the state is stable in a rubidium condensate regardless 
of the value of po. Fig. [5] illustrates how the stability 
depends on the value of q and the population pQ. We 
plot only the case = as it gives the fastest-growing 
instabilities and the smallest size of a stable condensate. 
We found numerically that the stability properties are 
independent of the value of 9. We have set 9 = m the 
calculations described here. If g > 0, the amplitude w' 
of the short-wavelength instabilities is suppressed as po 
increases. This can be understood with the help of the 
energy functional 



■2;0;-2 



:9infz +8qil- Pq). 



(38) 



li q > 0, the energy decreases as po increases. Therefore 
there is less energy available to be converted into the 
kinetic energy of the domain structure. From the top 
row of Fig. [5] it can be seen that Eq. ([57)) gives an upper 
bound for the size of a stable condensate also when the 
value of Po is larger than zero. If g < 0, the state is stable 
at Po = 0. The bottom row of Fig. [5] shows that now 
V'2;0;-2 bccomcs morc unstable as po grows. Increasing 



6 



po leads to an increase in the energy E2-o;-2- This can 
be converted into kinetic energy of the domains. Fig. [2] 
shows that Eq. ((37|) gives an upper bound for the size of 
a stable condensate also in this case. 



Nonzero tpi and 



As the next example we consider a state of the form 













(39) 



Also for this state (Fx) ~ (Fy) — and (Fz) = nf^ and 
therefore the Hamiltonian and time evolution operator 
are given by Eqs. ^ and respectively. The Bogoli- 
ubov matrix reads 



(40) 



Here Bi._i{t) time-dependent with period T = n/q and 
Bf._i is time-independent. The eigenvalues of Bf._^ are 



2,3,4 



± 



Efc (ek + (50 + 5i)n 



± n^/ {go - giY + 4goffi/^ 



1/2 



(41) 



Now ++, H — , — h, and correspond to wi, a;2, ^3, and 

a;4, respectively. These are all gapless modes. For ru- 
bidium the eigenvalues are real. In Appendix C we show 
that uji^2 are density modes and 0^3.4 are magnetization 
modes. 

We now turn to the eigenvalues of Bi._i. At q = 
Bi._i becomes time-independent and the eigenvalues are 



5,6 



±efc, 



± Sini/ (6efc + gmf^Y - 8ekf'^{iek - gm) 



(42) 



1/2 



(43) 



For rubidium these are all real. One of the eigenvalues 
/ia;4 5 has an energy gap \ginfz\. These eigenvalues de- 
scribe spin-magnetization modes. 

For non-zero q the stability can be analyzed using Flo- 
quct theory. As in the previous section, the fastest grow- 
ing instabilities are obtained at fz = 0. This case can 
be studied analytically by changing basis as described in 
Appendix [Cl The eigenvalues for the case /z = are 



^7,8,9,10 



-iq±^{ek+iq){ ek + 2gin + 3g), 



3g± 



(efc + qf + A{ekgin + q^) 



± A J [g2 + ek{gin + q)^ - Sgmq^el - q^) 



1/2 



(44) 



(45) 



These are gapped excitations with a magnetic field de- 
pendent gap. In more detail, at = we get 
^"^5.6 — ~3q ± ^iq{2gin + iq) and ^1^7,8,9,10 = 3(7 ± 



5g2 ± Aq^J q{3gin + q). For positive q, the fastest 

growing instability is determined by Wg and is located 
approximately at ek = max{0, — 3 -I- q}. For negative q 
there are three local maxima for a;'. The one with the 
largest amplitude is given by u}\ and w^q and is located at 
ek « max{0,q^(|g|-l)/(g^-|-|g|-|-l)}. The second-largest 
is given by and is at w max{0, 3|(?| -I- 1}. Finally, 
the instability with the smallest amplitude is related to 
uj^ and is at ek « max{0, 29(10|g| - 1)/100}. In Fig. © 
we plot the behaviour of w' for q = 6 and q ~ —3. From 
Eqs. and (^5]) it can be seen (see also Fig. ^ that 
the state is stable if the size of the condensate is smaller 
than 



Al;-1 



^2m{2\q\-q) 



(46) 




A. 



8 10 



FIG. 3. The positive imaginary part oj' of the eigenvalues 
(|44p and (|45p for g = 3 and q = —6. The unit of q and ek is 
\gi\n, and that of cj' is \gi\n/h. We have chosen = as it 
gives the fastest-growing instabilities and the smallest size of 
a stable system. The solid and dashed lines correspond to ojg 
and cjs, respectively, while the dotted line gives cjy and lu'iq 
[see Eqs. and ((45)) ]. 



C. Nonzero tp2 and i/)o 



As the final example we consider a state 

V'2;0 = Jt7 , 0<fz<2. 

V ^ 

V / 



(47) 



As for other states considered in this article, now (Fx) = 
(Fy) = and (Fz) = nfz and the Hamiltonian and time 
evolution operator are given by Eqs. ([8|) and ([9]), respec- 
tively. We note that the stability properties of the states 
■02;o and ?/'o;-2 = ^/r^(0,0, y/2 - fz,0, y%)/V2 are sim- 
ilar. Therefore the latter state will not be discussed in 
more detail. The Bogoliubov matrix of '!/'2:0 reads 



S2;0(t) 



-^2:0 ' 



' Bio ' 



BiM), 



(48) 
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where only Bj.q is time-dependent (with period T 
ir/q). The eigenvalues of and Bq.2 are 



fiwi 2 = ± Cfc, 



3,4,5,6 



± 



1/2 



(49) 



(50) 



These are gapless excitations. In Appendix |D] we show 
that uj'iA correspond to magnetization modes, while 
wi,2,5,6 are density modes. For rubidium, these are all 
stable modes. 

After a suitable change of basis the Bogoliubov matrix 
B|.q becomes time-independent, see Appendix (D] The 
eigenvalues of the new matrix are found to be 



1 



fiw7,8,9,io = '^~^\J *i =■= V (2efc + ginfz + 2gr)s2, (51) 



state 



V'2;-l, 



1p2;0;-2 



tp0;-2 



> 



< 



> 



< 



> 



< 



> 



< 



stable size 



27rft 



\/2m\q\ 



yj 6m I g 



\/6"i|9| 



fastest-growing instability [ek] 



2q — gin 



|(2T/.)|ffiln 



3g — gin 



3gin 



Q (k|-9l") 
"J^+klsin + Csi")^ 



-2-F0.9g + 0.04|/,|(l + q) 



lgl-3 + 1.3!M-0.16/; 



TABLE I. Summary of the results. Stable size gives the 
largest possible size of a stable homogeneous condensate 
and the fastest-growing instability indicates the approximate 
value of Efe corresponding to the fastest growing instability. If 
q is such that the ek given in the table is negative, the fastest- 
growing instability is at tk = 0. On the second line of the 
table, the — sign holds for V'2;-i and the + sign for i/'i;-2. 



where 

si = 2el + [gmf^f + 4efe[(3 - f,)gin + q] 

- 8fz9inq + 2q{6gin + 5q), 
S2 = fzigmfmek + q)- lOefc/, - ISqf, + gmf^] 

+ S2q^[q + efe + 3gin{2 - fz)] - IGginqekfz- 

Now ++, H — , — h, and are related to ^7, ojg, wg, and 

wiO; respectively. These are gapped excitations and cor- 
respond to spin-magnetization modes. These modes can 
be unstable for rubidium; an example of the behavior of 
the positive imaginary component of 1^7,8,9,10 is shown in 
Fig. HI 

An upper bound for the size of a stable condensate is 
the same as in the case of V'li-ij see Eq. With the 

help of Eq. (|5T|) it can be seen that the fastest-growing 
instability is approximately at = max{0, —2 -I- 0.9q + 
0.04/^(1 + q)} when q > and at Ck = max{0, |g| - 3 + 
1.3/^ - 0.16/2} when q<0. 




FIG. 4. The positive imaginary part cj' of the eigenvalues (jSlf) 
for q — 5 and q — —2. The unit of q and ek is \gi\n and that 
of uj' is \gi\n/h. 



V. CONCLUSIONS 

In this article, we have studied the dynamical stability 
of some nonstationary states of homogeneous F = 2 ru- 
bidium Bose-Einstein condensates. The states were cho- 
sen to be such that the spin vector remains parallel to the 
magnetic field throughout the time evolution, making it 
possible to study the stability analytically. The stabil- 
ity analysis was done using Bogoliubov approach in a 
frame of reference where the states were stationary. The 
states considered had two or three spin components pop- 
ulated simultaneously. These types of states were found 
to be stable in a rubidium condensate in the absence of 
a magnetic field, but a finite magnetic field makes them 
unstable. The wavelength and the growth rate of the in- 
stabilities depends on the strength of the magnetic field. 
The locations of the fastest-growing instabilities and the 
upper bounds for the sizes of stable condensates arc given 
in table m For positive q, the most unstable state, in the 
sense that its upper bound for the size of a stable con- 
densate is the smallest, is ip2;0;-2- However, this is the 
only state that is stable when q is negative. For q < 0, 
the states giving the smallest size of a stable condensate 
are '0i;-i and 4'2;0- 

In comparison with F = 1 condensates, the structure 
of the instabilities is much richer in an = 2 conden- 
sate. In an _F = 1 system, there is only one type of 
a state whose spin is parallel to the magnetic field. The 
excitations related to this state can be classified into spin 
and magnetization excitations (23j . In the present sys- 
tem, there are many types of states which are parallel to 
the magnetic field; we have discussed six of these. In ad- 
dition to the spin and magnetization excitations, there 
exist also modes which change spin and magnetization 
simultaneously. The increase in the complexity can be 
attributed to the number of components of the spin vec- 
tor. 
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The stability properties of the states discussed in this 
article can be studied experimentally straightforwardly. 
These states had two or three non-zero components, a 
situation which can be readily achieved by current exper- 
imental means [2^. Furthermore, the stability of these 
states does not depend on the relative phases of the pop- 
ulated components, making it unnecessary to prepare 
states with specific relative phases. 



Appendix A: Eigenvectors of B^.-i 

Here we give the (unnormalized) eigenvectors of i32 _i' 
i?!--!, and B2._i- Unlike -B|._i, the operators B2.^i and 
-Bl'-i depend on the magnetic field. The eigenvectors of 
the latter two arc given at g = 0. The operators i?2 _iJ 
-B2 _i, and i?2 _i not be given here explicitly as they 
can obtained straightforwardly from Eqs. and ([T^. 
However, we give the bases with respect to which these 
operators and their eigenvectors are defined. The matrix 
B2 -1 is given in the basis {bi, b4, bg, bg}, where b^ is a 
10-component vector with the j:th clement equal to one 
and all other elements equal to zero. The eigenvectors of 
B2-1 are 

= ((efc -t- fiwj)aj, {ck + hwj), {ck - fihjj)aj,ek - hwj), 

(Al) 

where j ~ 1, 2, 3, 4 and 



^ _ + fz 90 + 4ffi ^ Qek9o9in{2 - f^) 



' Y 2 - /, go - 2gi ' (.9o + 4gi)[e2 - (hu:,)^] 

(A2) 

The matrix Sf .i defined in the basis {b2, b5, bg}. At 
q = the eigenvectors are 



X5 =(V2-/.,-v/l + /.,0), 

X, =(.91 jV(2 -/.)(! + ./.), 51^^(2 - /.), 



(A3) 



{ek+29in- huj)), .7 = 6,7. (A4) 

By defining -B|._i with respect to the basis 

{b7,bio,b3} we get B2'_i = ^{B2;-i)*- Therefore, the 

eigenvectors of Sf -^i can be obtained from those of -Bf -^i 
by complex conjugating. 

Appendix B: Eigenvectors of i?2;0;-2 



m 



r6 

^2;0;-2 



The operator B2 Q -2 appearing in Eq. ([33|) is given 
the basis {bi, ba, bs, bg, bg, bio}. The eigenvectors of 
corresponding to cji.2 are 



xi =(VPoP-2, -e*^V/'2P-2, 0, 0, 0), 

X2 =(0, 0, 0, VpoP^, -e*^ VP2P^, VP0P2). (BI) 



These are magnetization modes as they change the mag- 
netization but not the spin direction. The exact eigen- 
vectors corresponding to a;3^4.5.g arc too complicated to 
be given here. Therefore we approximate 91 w (for 
rubidium 91/90 « 0.01) and obtain 

X3,4 ~(VP2' e'^VPo, y^P^, 

-7±\/p2, -7±e~*''VPo,-7±VF^)' (B2) 
X5 «((2 - DVpJ' -e^'/.VA^, -(2 + /.)Vp^,0,0,0), 

(B3) 

xg «(0, 0, 0, (2 - /,)VpI, -e^'f,^, -(2 + /.) VP^), 

(B4) 



where 



P±2 =i(2-2po±/.), 



7± 



1 

gon 



?fe + gon ± \J tk (e/c -f 2go»^) 



(B5) 
(B6) 



Of these X3.4 arc density modes and X5 g arc magnetiza- 
tion modes. 

The operator ^^20-2 given in the basis 

{b2, b4, by, bg}. B\.q._2 is time-dependent, but at 
Po = the time evolution determined by B\.q._2 can be 
solved analytically. With the help of the unitary basis 
transformation 



(B7) 







g-3it(j 





\ 


1 








g-3itg 




71 










_g3itij 











/ 



wc obtain a new Bogoliubov operator 



r4 

^2:0;-2 



= Ws4 v + ih{-V^]v, (B8) 



which is time-independent. The eigenvectors of this op- 
erator are 

X7.8 =(2/1^4,5 - .giTi/z - 6g, 



2e, + [2+ VWj]ffin-6g,0,0), (B9) 
X9,io =(0, 0, 2huj(;j + ginfz - 6g, 

2e, + [2+ ^4-./|]gin-6<7). (BIO) 

These modes change both magnetization and spin direc- 
tion. 



Appendix C: Eigenvectors of Bi--i 

The operator Bf._i is defined in the basis 
{b2, b4, by, bg}. Thc eigenvectors arc (in the limit 
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.91 = 0) 



Appendix D: Eigenvectors of B2-o 



xi = (v/TT/T, VW^, -7- Vi+TT, -7- VW.)> 

(Cl) 



X2 = (-7- V^TTz, -7- v^wT, yTT/:, v/wT), 

(C2) 
(C3) 
(C4) 



X3 = (v/W^,-Vi+X,o,o), 

x4 = (o,o,v/wr,-yTT7:). 



Here 7± is defined as in Eq. (|B6[) . Clearly xi^2 are density 
modes and X3^4 are magnetization modes. 

The operator -B2;0;-2 is defined in the basis 
{bi, ba, bs, bg, bg, bio}. Here we give the eigenvectors 
at /, = 0. 



The operators i?|.o,-B|.o, and i?|.Q are defined in the 
bases {bs, bio},{bi, bs, bg, bs}, and {b2, b4, by, bg}, re- 
spectively. The eigenvectors of B^.g and B|.q read 

xi,2 =b5,io, (Dl) 

X3,4=("7±yA,-7±VWr,//^,v/2^), (D2) 

X5=(v/2~7:, 777,0,0), (D3) 

X6=(0,0,v/2-/.,v^). (D4) 

Here "f± is defined as in Eq. (|B6p and the index 3 cor- 
responds to 7+ and the index 4 to 7_. Of these X3.4 
arc density modes and xi.2,5.6 are magnetization modes. 
When calculating X3.4^5.g we have set 171 = 0. 

The time-dependence of the operator i?|.g can be elim- 
inated with the help of the basis transformation 



X5,6 = (ffc +gin- fia;5,6, 51*^7 0,0,0,0), (C5) 
X, = (0, 0, a„ /3j, 7,, <5,), J = 7, 8, 9, 10, (C6) 



where aj , /3j , , Sj are too complex to be given here. 
These modes change both spin direction and magneti- 
zation. 










e-**9 


\ 


1 




























1 







; 



(D5) 



The eigenvectors xy^g.g.io of the resulting time- 
independent operator describe spin-density modes and 
are too complicated to be given here. 
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